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We calculate, by means of fiuctuational electrodynamics, the thermal emission of an aperture sep¬ 
arating from the outside, vacuum or a material at temperature T. We show that thermal emission 
is very different whether the aperture size is large or small compared to the thermal wavelength. 
Subwavelength apertures separating vacuum from the outside have their thermal emission strongly 
decreased compared to classical blackbodies which have an aperture much larger than the wave¬ 
length. A simple expression of their emissivity can be calculated and their total emissive power 
scales as T® instead of T"* for large apertures. Thermal emission of disk of materials with a size 
comparable to the wavelength is also discussed. It is shown in particular that emissivity of such a 
disk is increased when the material can support surface waves such as phonon polaritons. 

PACS numbers: 44.40.-|-a,05.40-a,78.67.Pt 


INTRODUCTION 


Since the end of the 19th century and the work of Max Planck, it has been known that thermal emission of radiation 
follows universal laws. For instance, the emissive power of a body at temperature T cannot exceed the value given by 
the so-called Stefan law, that reads as i/°(T) = aT^, with cr = 5.67 10“® W.m“^.K“'*. Another feature is that the 
thermal emission spectrum is broadband and peaked around Am (given by the Wien law AmT = 2898/im.K), with 
a typical bandwidth of a few Am- However, theoretical models based on a fluctuationnal electrodynamics formalism 
have shown that thermal emission could deviate from the above mentioned behaviors when the length scales involved 
are small compared to the typical wavelength of the emitted radiation. For example, when two heated bodies are 
separated by a small gap, radiative heat transfer surpasses that predicted by classical formulas, due to the coupling of 
evanescent modes on the surface of each body [l|,[3- Heat transfer is enhanced in this case, and can even be dominated 
by transfer through modes at specific frequencies, especially when the materials exhibit resonances such as surface 
phonon or surface plasmon polaritons Moreover, micro or nanostructured surfaces, such as periodic gratings, can 

scatter the thermaly excited evanescent waves into the far field, which substantially changes the emission properties. 
This mechanism has paved the way towards the design and fabrication of coherent thermal sources exhibiting both 
temporal and spatial coherence Q • Another way to couple the near held and the far held is to use the tip of a Scanning 
Near-Field Optical Microscopy and bring it at a submicron distance from the heated surface. The thermally populated 
evanescent modes can be coupled to a detector in the far held by scattering at the tip. This process underlies the 
principle of Thermal Radiation Scanning Tunneling Microscopy 1413 (TRSTM), an imaging technique among others 
[ll| that uses thermal radiation to perform imaging and spectroscopy of subwavelength structures. 

The purpose of this paper is to explore another aspect of thermal emission at subwavelength scale. We study the 
conceptually simple situation of thermal emission by an aperture. Note that the problem could be addressed using 
the reciprocity theorem of electromagnetism. Indeed, in the theory of thermal radiation, it is known that reciprocity 
is the foundation of Kirchhoff’s law, stating that the emissivity of a material equals its absorptivity. This means that 
knowing the absorption efficiency (5abs(w) of a body at a given temperature T, the thermally emitted flux by this 
body at the same temperature and at frequency oj is given by 121 


= Qahs(w) 


huj^ 


4TT^c^[exp{huj/kiiT) — I] 


( 1 ) 


where h is the reduced Planck constant and kj, is Boltzmann’s constant. Therefore, the knowledge of the light 
absorption properties of an object at a given frequency allows one to deduce its thermal emission properties. For 
example, a sphere of a homogeneous material will emit according to Fq. o with Qabs given by the Mie theory 0. 

In this paper, we address the problem from a different point of view. We use fluctuationnal electrodynamics in 
order to compute directly the thermal emission by an aperture. The principle of the approach is the following. In a 
body at local thermal equilibrium, temperature initiates fluctuating currents that radiates an electromagnetic held 
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FIG. 1: Geometry of the model system. From the electromagnetic held in the plane z = zq, one deduces the held and the 
radiated power in any plane at a distance z. 


14j . Thermal currents are characterised statistically by a correlation function given by the fluctuation-dissipation 


theorem. Radiation by these currents is calculated by solving Maxwell’s equations in the specific geometry, as in a 
standard antenna radiation problem. Note that in the specific case of an aperture, the emitted heat flux is given by 
the flux of the Poynting vector through a plane parallel to the aperture, allowing us to connect this flux to the Wigner 
transform of the electric field spatial correlation function 15l-ll7j. These spatial correlations are directly computed 
in fluctuationnal electrodynamics [l^, [I^ . We first focus on the simple case of an aperture separating vacuum at 
thermal equilibrium from the outside. Then the formalism is also applied to the case of an aperture separating a 
material supporting resonant surface waves at thermal equilibrium from the outside 


EMISSIVITY OF AN APERTURE 

The system considered here is depicted in Fig. [TJ A semi-infinite material at temperature T fills the half-space 
z < 0, on top of which a mask with a transmission function t(R) is placed in a plane z = zq, where zo —>■ 0. Thermal 
radiation is emitted by the material through the mask, and the radiated power is calculated in a plane at a distance 
z through the evaluation of the flux of the Poynting vector across this plane. 

For monochromatic fields, the complex amplitude of the electric field E(r) in the plane z can be written as a 
plane-wave expansion in the form 


E(r) = / (2) 

J 47r^ 

where k = (K, 7 ) = Ku± + , r = (R, z) and 7 ^ -I- = /cq with kg = ojIc = 2tt/X. The amplitude E(K, zq) of 

the plane waves in this expansion is the Fourier transform of the field in the plane z = zq, and reads as 

E(K,zo) = J 'E{Ro,zo)e-^'^-^°d^Ro ■ (3) 

The power 4>{uj) radiated in the far field is defined as the flux of the Poynting vector through the plane z. For 
monochromatic fields, the time-averaged Poynting vector is S(r) = l/25i [E(r) x H*(r)], where H(r) is the complex 
amplitude of the magnetic field and the superscript + stands for complex conjugate. Using the Maxwell equation 
V X E = jw/ioH and the plane-wave expansion of the electric field, one obtains 




2fiQUJ 




7|E(K,zo)|" 


cPK 

dTT^ 


( 4 ) 


Note that the integration is restricted to propagating waves, i.e, waves for which K < kg since 5 R( 7 ) = 0 when 
K > kg. In this case, this integration can also be understood as an angular integration on the upper hemisphere of 
the wavevector k with constant modulus |k| = uj/c. 
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Equation dH shows that the knowledge of the field in the plane z = zq permits an explicit calculation of the 
radiative flux emitted in the far field. In our model, this field can be understood as the field radiated by the semi¬ 
infinite medium and transmitted through the aperture. Denoting by the field right before the plane of the 
aperture, and describing the aperture (or actually any scattering object placed in the plane z = Zq) hy a. transmission 
matrix Tij(K, K'), one can write the field in the plane z = zq as 

E,{K,Zo) = I r,,(K,K')£;r(K',2o)^ • (5) 


Inserting this expression into Eq. @ leads to 


= 


1 


32/ioW7r® 




(K, K')t4(K, K")K”“(K', zo)E“=*(K", Zo)d^Kd^K'd^K'‘ 


( 6 ) 


The incident field can be calculated as the field radiated by the semi-infinite material in absence of the aperture 
(this is the simplest model, a self-consistent calculation being outside the scope of the present study). This field is 
linearly related to the thermally excited electric currents inside the material, through a relationship of the form 


"(r) =ifioU} / d^r'Gim{r,r')jm{r') 


(7) 


where G is the tensor Green function that describes the electrodynamic response of the semi-infinite material and j 
is the electric current density. The Green function in this geometry can be written as a plane-wave expansion that 
involves the Fresnel transmission factors at the interface z = 0 between the medium and vacuum [^ . According to 
this expansion, the incident electric field reads 


Er(r) 


Stt^ 


d^r'd^K 

72 


[eiSt 2 iS jm(r') 


I Er(K,zo) 


d^K 

47r^ 


( 8 ) 

(9) 


where s = K/|Ar| xe^, pf = [K'^e^ — yiK^ex — liKy^y\ /(nikoK), and and are the Fresnel transmission factors 
for s and p polarization, respectively (^ . By identification, one obtains the expression of the Fourier transform of 
the incident field in the plane z = Zq : 


Er(K,zo) = 


-poUJ 


d^r' r „ g ^ 

- H" 

72 




( 10 ) 


The thermally excited currents are fluctuating fields, that are describes statistically. In order to compute fluxes, one 
needs second order quantities. The spatial correlation function of the currents in the material at thermal equilibrium 
is given by the fluctuation-dissipation theorem 




eoS[e(aj)]aje(u;,T) ^^^^^^ 

TT 


r')S{uj-cj') 


( 11 ) 


where the brackets denote an average over thermal fluctuations, 0(a;,T) = huj/[eyip{huj/kbT) — 1], and e(a;) is the 
dielectric function of the medium. From Eqs. ©, (nni) and dm, one obtains the following expression of the thermally 
radiated flux 






d^Kd^K'y{K) 


jft(72(K')) 

|72(K')P 


(K,K')r4(K,K') 




where 


M^(K) 


\3£ 


Kl -K,Ky 0 \ 
-KxKy Kl 0 
0 0 0 / 


and 


MP(K) 


|,p ,2 \12?+K^ 

' |n2P|ni|2fc4if2 


\y\^Kl \y\^Kx,Ky \ 

\y\^K,Ky \y\^Kl -yKyK^ 

-y*K,K^ -l*KyK^ ) 


( 12 ) 


(13) 


(14) 
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Note that this final expression is restricted to positive frequencies only (as is usual in radiative transfer), which 
implicitely assumes that all fields in the derivation have been replaced by their analytic signals (in practice this 
results in an extra factor of 4, see @ for details). 

This expression of the radiated power appeals for the definition of an effective emissivitty. Indeed, in the framework 
of geometrical optics, the emitted flux by an object with surface S is usually written in the form 


</.(cu,T) 


Q{uj,T)uj^ 


(15) 


where e is by definition the emissivity of the object. From Eq. (1121) one can dehne the effective emissivity of the 
aperture [or of any scattering object defined by a transmission matrix Tij (K,K')] as 


^eff — 


Stt^/cqS' , 




, ^[72(K')] 

|72(K')P 




(K,K')r4(K,K') M;,(K')+Mj;(K') (16) 


This is the general expression of the emissivity of an aperture defined by its transmission matrix Ty (K, K'). It involves 
a double integral the transmission matrix over all parallel wavevector. Integration over K is limited to propagative 
waves such as K < fcp, whereas integration over K' includes a priori both propagating {K' < kg) and evanescent 
waves {K' > kg). The contribution of evanescent waves to the radiated flux in the far field results from a scattering 
process. The thermally excited evanescent waves with large wavevectors K' > fcg a-re scattered into propagating waves 
with K < k^hy scattering at the aperture. Another feature of the expression of the effective emissivity is that the 
material and geometrical resonances are contained in the integral both in the transmission matrix Tij(K, K') and the 
Fresnel transmission factors. Finally, note that due to reciprocity, the expression of the emissivity can also be seen as 
that of the absorption cross-section normalized by the geometrical cross-section S. 


APERTURE IN VACUUM 


As the simplest example, we consider the case of blackbody radiation in a vacuum at temperature T transmitted 
through an aperture in an opaque screen. In the general model derived in the preceding section, this amounts to 
considering a material with transmission factors t'* and equal to unity. The radiative heat flux coming out from 
the aperture can be calculated analytically in two asymptotic cases. The first case corresponds to an aperture with 
a radius cq much larger than the typical thermal wavelength. Under this assumption, one can make use of the 
Kirchhoff approximation in which the field equals the incident field in the aperture and vanishes outside. The limit 
of validity of the Kirchhoff approximation is estimated to be k^r^ ~ 6, which corresponds to an aperture radius 
on the order of the wavelength 2l|. Under this assumption, the transmission matrix is reduced to a scalar so that 


, (K, K') = 5^JT{K - K'), where 


T(K) = / r(R). 


zK.R j 2 




(17) 


and T(R) = 1 inside a circle of radius vq (the aperture) and r(R) = 0 outside. An explicit calculation leads to 
r(K) = ^ dip Jo° 

Inserting this expression of the transmission matrix into Eq. (1121) allows in principle to calculate the radiated flux. It 
is however easier to rewrite the flux as 


<(>(w,r) = 


0 (w,r) ^(K) 


d^Kd^K'd^Ti'd'^R" 


IGtt® J 7(K' 

and to perform the change of variables m = (R' 4- R")/2 and d = R' — R", leading to 

e{oj,T) f yiK) 


</>(u;,T) = 


^T(m -f d/2)r(m - d/2)e-^^^-^''> ‘^d^Kd^K'd^md^d . 


IGtt^ J 7(K') 

Since the product T(m -|- d/2)T(m — d/2) is independent on the variable m, the integration over m gives 

d 


J r(m -I- d/2)T(m — d/2)dm = 7rrQtU(d) = ttcq x 


arccos ■ 


2rn 


d 

2ro 


' 1 - 


A 

2ro 


(19) 


( 20 ) 


( 21 ) 
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Using spherical coordinates with angles 9 and ip, one can write k = (K, 7 ) = fco(sindcos(/3,sin0sin(/5,cosd) and 
transform the integral into 




e(a;,r) 

IGtt® 


Trro fcp 


cos^ 0 sin 6»e ® sin 0' ®V (d) dOdipde' dip'd^ d 


( 22 ) 


which, after integration over azimuthal angles, gives 




e(cu,T) 

IGtt® 




J cos'^ 9 sin OJo^kod sin 9) sm9'Jo{kQd sin 9') W{d)d9d9'd‘^d . 


(23) 


Integration over 9 and 9', knowing that d extends over a disk of radius 2ro, leads to 


0 ( 0 ,, T) = 




2 r2koro 

7rrQ2 / IU(u/fcQ) sinu-F'(M)du 


/o 


p2koro 

= 0 *^( 0 ;) / W{u/kQ)2 sin uF{u)dv 

Jo 

= </'°(w)eS{{(w) 


(24) 


where F{u) = (sinu — ucosu)/vP‘. The last expression defines the effective emissivity sl^{uj) at frequency u; of a 
blackbody of circular radius tq. 

When the aperture is large compared to the wavelength, thermal emission corresponds to a blackbody. However, our 
result shows that the emissivity of an aperture is smaller than 1 if the aperture size is on the order of the wavelength. 
Pushing the Kirchhoff approximation at its limit kor^ = 6, we obtain ~ 0.84. This can be easily understood since 
waves with wavelengths on the order or smaller than the aperture size can hardly be transmitted. The aperture acts 
as a high pass filter, reducing the contribution of low frequency waves, which is a feature of the underlying diffraction 
process. However, it is known that the Kirchhoff approximation breaks down when the aperture size becomes smaller 
than the wavelength (2ll42^ . typically when koro < 6. Bethe [ 2 ^ and Bouwkamp [ 2 ^ have indeed shown that the 
transmission through a small hole is actually weaker than that predicted by the Kirchhoff approximation. The problem 
addressed by Bethe and Bouwkamp’s theory is that of transmission through a hole in a perfectly conducting screen. 
By introducing fictitious magnetic charges and currents in the diffracting hole satisfying boundary conditions on the 
screen, their theory allows one to calculate the scattering cross-section and the transmission matrix r(K,K') in the 
regime koro <C 1. One ends up with 


‘ik 


2 _ 64 2 6 cos^ 


'K,K')\^ = -k, 


0'0 


cos^ 9 


(1 — sin^ dcos^ ip) 


(25) 


for s polarization, and with 


' ik 


, 2 ^^ 7.2 qCOs‘^ 6sin^ 6{cos‘^ ifl/4:COs‘^ if') — sinO cosifsinO' 

(->^5 j| — ~::rknTn 


^ 0 ' 0 " 


cos^ 9 


(26) 


for p polarization. Let us note that the transmission matrix is here limited to propagative waves {K,K' < kg). 
Inserting these two expressions into Eq. 1121 one can perform the integration over incoming and outgoing wavevectors 
(K an K'), which for propagating waves amounts to integrating over 9, ip, 9' and (p'. This leads to the following 
expression of the radiative thermal flux emitted by a subwavelength hole : 




16 fcgrg 
27 7r3 


0 (cu,T) 


64fcgrg 

277r2 


4P{u:,T) 


£e//(/>°(w,T) . 


(27) 


This result shows that the effective emissivity of a subwavelength hole is Seff = 64(fcgrg)^/(277r^). As expected, this 
emissivity is smaller than that predicted by the Kirchhoff approximation, which predicts a scaling in fcgTg. Note that 
the scaling in fcgr-g that is obtained for a subwavelength hole is consistent with that expected for Rayleigh scattering 
(ie. scattering by particles much smaller than the wavelength). This result confirmes that small apertures behave as 
high-pass filters regarding thermal emission. 

Expression (lT7l) gives the radiative flux at a given frequency w. If the condition /coro 1 is satisfied on the full 
spectral range covered by thermal emission (typically Am/2 < A < 5Am in terms of wavelengths), the spectrally 
integrated flux can be calculated, and reads 


(f) = 



16 fcgrg 

27 7r3 


0(w, T)duj 


405c^F 


X Trrg . 


(28) 
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It is interesting to note that instead of following the usual law of free-space blackbody radiation, the power emitted 
by a subwavelength blackbody follows a T® law. This means that for a given aperture size tq , when the temperature is 
decreased so that is larger than tq, the thermally emitted power decreases drastically, much faster than predicted 
by the usual Stefan-Bolztmann law. For example, a hole with tq = l//m at 77 K (liquid Nitrogen temperature) has an 
emissive power of 1.99 W.m“^ according to Stefan-Boltzmann law, and of 4.75x 10“"^ W.m“^ according to the law 
derived in this paper using the Bethe-Bouwkamp theory. Finally, let us remark that deriving an analytical expression 
of the emissivity in the intermediate regime fco^o ~ 1 is out of reach. In that case, one should follow approaches 
that have been used, for example, to address the problem of extraordinary transmission through subwavelength holes 
and compute the absorption efficiency, that directly leads to the emissivity according to Kirchhoff’s law. 


APERTURE EILLED WITH A MATERIAL 


In this section we address the thermal emission by an aperture when the medium occupying the half-space z < 0 is 
a real material (see the geometry in Fig. 1). This problem cannot be solved in its full generality since there is no exact 
expression of the transmission matrix r valid for any material. However, the Kirchhoff approximation can be used as 
long as fcofo ~ 6, and we limit the study to that regime. This will allows us to highlight interesting phenomena that 
occur when the aperture size approaches the wavelength. Under the Kirchhoff approximation, the emitted radiative 
flux reads 

(j){L0) 


X 


^2/coro 

A°(w) / W{u/k{))uF(u)du 
Jo 


KJo{Ku)dK 
3 a/I — 


(2 - |r-|2 _ |^P|2) + r 2Ac Jo(At»)d 3t ^ 

Jl 


(29) 


where k = K/k^. This expression contains two contributions: the propagating wave contribution for k < 1 and the 
evanescent wave contribution for k > 1. Let us first check that from Eq. (l29ll one recovers the classical expression of 
the radiative flux when the aperture size is much larger than the wavelength. For a circular aperture, this corresponds 
to the condition kovo 1. Let us note that W{u/ko) decreases smoothly from 1 to 0 when u/ko varies from 0 to 2ro- 
F{u) decreases fastly to 0 when u is large compared to 1. When /cpro ^ there is a domain in which u ^ 1 and 
u <g; 2fcoro- In this domain, the upper bound of integration over u in Eq. ((29ll can be replaced by oo, and W{u/kQ) 


can be replaced by 1. Noting that uF{u)Jo{Ku)du vanishes if k > 1 and Vl — if /t < 1 27|, one retrieves that 


there is no contribution of the evanescent waves to the emitted flux for large apertures. Moreover, the expression of 
the emitted flux equal the classical expression 


(p = (j^cias = [ KdK{2 - - |r^p) (30) 

where the integral represents the emissivity of the material. Note that this emissivity is equal to 1 when the Fresnel 
reflection factors vanish, i.e. in the vacuum blackbody radiation limit. 

In the regime where fepro is not large compared to one, the contribution of the evanscent waves is no more negligible, 
and one has to integrate Eq. (1291) numerically. An interesting situation is that of a material supporting surface waves, 
such as SiC, at the limit of validity of the Kirchhoff approximation in terms of aperture size. In Fig. [21 the effective 
emissivity {i.e. 4>/4P) is plotted versus frequency around the surface-phonon polariton resonance of SiC which occurs 
for A = 10.6 /xm. For an aperture with radius rp = 100 /xm filled with SiC, the emissivity is the same as that obtained 
for a massive material. It is close to one in a broad specral range, except close to the surface-polariton resonance 
for which the material is very reflective. For a radius rp = 10 ^m, the emissivity is enhanced in the spectral domain 
where SiC supports surface polaritons. These surface polaritons are thermally excited and scattered by the aperture, 
which adds new channels for far-held thermal radiation. One can even observe an effective emissivity larger than one 
around the surface polariton resonance frequency. This means that the thermal emission of the aperture is larger that 
the blackbody emissive power multiplied by the geometrical cross-section. A radiometric interpretation is that the 
effective aperture emission size is larger than it geometrical size. Using reciprocity (or Kirchhoff’s law), one can also 
understand that the emissivity is equivalent to an absorption cross-section, normalized by the geometrical section. 
It is actually well-known in scattering theory that scattering by nano-objects or nano-antennas such as nano-spheres 
or nano-cylinders leads to cross-section larger than the geometrical size. This is the so called antenna effect. Note 
however that when the surface considered for thermal emission becomes larger than the wavelength, there is no way 
that this emission can surpass blackbody emissive power. For example, it is not possible to make a macroscopic 
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FIG. 2: Emissivity vs angular frequency for a circular aperture filled of SiC with a radius of 100 um (plain) and with a radius 
of 10 /rm (Dashed). 


surface made of small aperture that overall would surpass blackbody limit. There is therefore no violation of the 
blackbody limit for macroscopic surface containing or not sub wavelength objects. 


CONCLUSION 

We have shown that thermal emission by a material can be substantially modified by confining this material to areas 
on the order or smaller that the typical emission wavelength. The confinement acts as a high pass filter, that changes 
the spectrum of thermal emission, as well as the value of the effective emissivity. In the case of a subwalength hole, the 
effective emissivity has been calculated using the Bethe-Bouwkamp model. It has been shown that in this limit, the 
emissivity scales as fco'''oi total emitted flux scales as T®, instead of the usual blackbody T'^ law. In the case 

of an aperture separating a material supporting surface modes (such as surface-phonon polaritons) from the outside, 
a contribution from evanescent waves scattered by the aperture generates an enhancement of the emissivity around 
the resonant frequency. From a thermal engineeering point of view, this study shows that the design of subwavelength 
scattering structures (the aperture being a simple example) could allow one to produce thermal sources with high 
spatial confinement and large efficiency at specific frequencies. The design of more complex structures would require 
an improvement of the theory to solve the full electrodynamic problem without requiring simple geometries or crude 
approximations. This could be done using numerical approaches already in use in nanophotonics, and in fluctuating 
electrodynamics such as discrete dipole approximation (DDA), finite-domain time difference (FDTD), or rigorous 
coupled wave algorithm (RCWA), to cite a few. 
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